Incommensurate DDW order 
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We consider various incommensurate (IC) order parameters for electrons on a square lattice which 
reduce to d^^z -density wave (DDW) order when the ordering wavevector Q — > (tt, tt). We describe 
the associated charge and current distributions and their experimental signatures. Such orders can 
arise at the mean-field level in extended Hubbard models. We compare the phase diagrams of 
these models with experiments in the underdoped cuprates, where (1) DDW order is a possible 
explanation of the pseudogap, and (2) there are experimental indications of incommensurability. 
We find various types of IC DDW and discuss their possible relevance to the physics of the cuprates. 
Our main finding is that IC DDW order is generally accompanied by superconducting order, but 
the magnitude of the IC wavevector can be small. A comparison with the analogous AF-ICSDW 
transition is given. 

PACS numbers: 



I. INTRODUCTION 

When a commensurate one-dimensional charge den- 
sity wave (CDW) is doped, the resulting charges may 
be viewed as defects in the CDW. If the interactions 
between them are sufficiently strong compared to their 
effective kinetic energy, then they will form an ordered 
lattice (which can be stabilized by a crossover to 3D or- 
der if there are many ID chains coupled together) of 
defects and the CDW will become incommensurate with 
the lattice. This can be more energetically favorable than 
simply doping holes into the rigid band formed in the 
presence of CDW order at a fixed wavevector. The rea- 
son is that, in the latter case, we must pay the single- 
particle energy gap of the CDW. In a two-dimensional 
-density wave (DDW) state, however, there are 
nodes in the order parameter, so it is relatively pain- 
less for doped holes to simply go to the nodes, which 
expand into Fermi pockets (in quasi-lD systems, how- 
ever, this might not be possible, so doping must lead to 
incommensurability 1 ). However, this cannot continue in- 
definitely if the DDW state is stabilized by approximate 
nesting, since the Fermi surface (FS) would eventually 
move away from the nesting wavevector. Thus, incom- 
mensurate DDW order is a strong possibility, at least over 
some range of hole dopings. In this paper, we explore this 
possibility. 

Near half-filling, the Hubbard model, the t — J model, 
and generalizations of these appear to have many phases 
which are nearby in energji^. Thus, small perturbations 
can strongly influence the competition between them. 
Consequently, it should not be too surprising that exper- 
iments on the cuprate superconductors have also uncov- 
ered evidence for a cornucopia of phases, particularly on 
the underdoped side of the phase diagram, which appear 
in particular materials and for certain values of the dop- 
ing level, temperature, magnetic field, etc£. Depending 
on the interfacial energies between these phases, one way 
in which their competition can be resolved is through 
the formation of stripes^ or other inhomogeneous pat- 



terns of microscale phase-separation. Neutron scattering 
and scanning tunneling microscopy experiments point to- 
wards this possibility^-' 6 . Furthermore, as the cuprates 
are doped, their Fermi surfaces evolve away from the 
commensurate nesting wavevector (it, 7r), so we would 
expect translational symmetry-breaking order parame- 
ters to occur at incommensurate wavevectors. Hence, it 
would seem important for any theory of the pseudogap 
to incorporate the tendency towards incommcnuration, 
which we take as further impetus to study incommensu- 
rate (IC) order parameters related to DDW order. 

In this paper we consider various ways in which d x 2_ y 2- 
density wave (DDW) order can go incommensurate. Be- 
fore turning to this discussion, we review the reasons for 
expecting DDW order to occur in the pseudogap regime 
of the cuprates. It has been suggested^ that competi- 
tion between DDW order and c£-wave superconductivity 
(DSC) can explain many of the anomalous properties of 
the pseudogap^. The DDW order parameter is given by 



=L c k+Q«x> = *$/(k) 



(1) 



with /(k) = cos(fc a ) — cos(k y ) and a commensurate (C) 
ordering wavevector Q = (tt, it). The particle-hole con- 
densate in Q breaks translation by one lattice spacing, 
parity, and time-reversal but is invariant under the com- 
bination of any two of the above symmetries. In real 
space it is represented by currents which alternate from 
one plaquette to the next. 

One motivation for considering DDW order is to 
explain^ the abrupt depletion of the superfluid density 
below the critical x c « 0.19 in Bi2212 and similar dopings 
in other materialsAi*^. Because the order has a e?-wave 
momentum dependence, its onset would explain ARPES 
measurements above T c in the pseudogapA?. Further- 
more, since DDW is a spin-singlet condensate, it would 
also explain the suppression of the spin susceptibility in 
NMR measurements and opening of a spin gap in in- 
elastic neutron scattering^. Also, the depletion of the 
c-axis conductivity 14 , whose contribution mostly comes 
from the antinodal regions, would also be naturally ex- 
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plained by the development of DDW order, which gaps 
these regions. Direct attempts to measure DDW order 
through neutron scattering have neither ruled it ouli^ 
nor unambiguously verified its presence, but there are 
intriguing suggestions that it may be present^. Thus, 
it is natural to ask whether DDW order is compatible 
with incommensurability. (For a discussion of the exper- 
imental signatures of DDW order, the reader is referred 
to the literaturo 7 i 10 i 13 i 17 i 18 i 19 i 20 i 21 . The case of DDW in 
La2- x Ba x CuOi is discussed in 2 ^. For numerical evi- 
dence of the existence of DDW in other systems, one 
is referred to^ 3 ..) 

The plan of our paper is as follows. In section ITT1 we 
develop some preliminary intuition about the energetics 
of IC density wave order by analyzing the DDW suscepti- 
bility of a Fermi liquid (both the nested and non-nested 
cases). In section ITTT1 we review some of the theory of 
commensurate (C) DDW order and show various ways 
to make it incommensurate. We also suggest the kind of 
microscopic physics that could lead to C and IC DDW, 
and in addition discuss possible experimental signatures 
of IC DDW. Finally, in section llVl we compute the mean- 
field phase diagrams of an extended Hubbard model. Our 
main results are that (a) IC DDW develops for a wide 
range of parameters, but (b) the IC wavevectors can be 
very small. Their relevance to the cuprates may, at most, 
be for a rather narrow set of dopings. We also offer a 
physical explanation of our results, which sheds some 
light on the nature of C-IC DDW transition in extended 
Hubbard models. We note that similar incommensurate 
order has been found in refs. 24 ' 25 . 



II. SUSCEPTIBILITY NEAR HALF-FILLING 

The easiest way to see that IC DDW phases could arise 
in effective models (which, in this paper, will, in turn, 
emerge from microscopic extended Hubbard models with 
correlated hopping terms2&; see section All CI for details) 
is to consider the bare susceptibility as a function of the 
incommensurability wavevector q, where Q = (ir, ir) + q: 



close to half-filling, where perfect nesting at half-filling 
was asumed. In fact, his results can be recovered by let- 



Xo(q) = ^ / (/(k)) 



2 n F (e k+ Q - m) - n F (e k - y) 



(2) 

where /(k) is defined after eq. ^ Here, for simplic- 
ity, DDW is considered without the presence of any 
other order. In eq. [3 np is the Fermi occupation 
number, and the energy dispersion contains nearest, t, 
and next-nearest neighbor, t', hopping parameters, e k = 
— 2t(cos k x + cos k y ) — 4t' cos k x cos k y . The DDW equiv- 
alent of the Stoner criterion is satisfied when the right- 
hand-side of this equation approaches the inverse of the 
pertinent coupling 1/g. This observation allows us to in- 
terpret the positions of the peaks of Xo(q) a s the wavevec- 
tors at which DDW order is likely to occur as we lower the 
temperature. A similar analysis was done by SchulaS 7 . for 
the case of AF and IC SDW order in the Hubbard model 



X 



120 



100 - 



DDW 



t =1 , t'=0, p=ioo/t 




FIG. 1: Susceptibility to DDW order. Away from half-filling, 
the peak of Xoil) is a t nonzero q. Here, q is measured along 
the zone diagonal in units of (n/a,n/a), where a is the lattice 
spacing. 

ting /(k) -y 1 in @ if we set f = 0. 

The support of (/(k)) 2 in the integrand in is at 
(0, 7r) and symmetry-related antinodal points of the FS 
since this is where the DDW order parameter is large. 
However, since these points are also the location of van 
Hove singularities, the regions around them also domi- 
nate the AF susceptibility. Thus, at least in the t' = 
case, one expects that whenever an IC SDW instability 
occurs, an IC instability for DDW should also occur. In- 
deed, for both orders, the peak of Xo(q) occurs at q ^ 
for low enough temperatures and sufficiently away from 
half-filling, e.g. if > 0(1), as shown in Fig. □ Note 
that the same analysis repeated with any other density 
wave or spin-density wave order would yield an IC insta- 
bility as long as the dominant channels of the modulation 
coincide with Fermi surface points which give large con- 
tributions to the susceptibility. For instance, ap^-density 
wave 8 - with wavevector Q = (0, n) + q has a peak in xa 
for q 7^ away from half filling at T ~ 0, in the case 
t' = 0. 

A natural question to ask is: to what extent does the 
IC DDW instability mimic the IC SDW one? To an- 
swer this, one needs to consider specific models as we 
do in IIII CI However, two differences can be anticipated 
even beforehand. First, the DDW order parameter has 
nodes, while the AF order parameter does not. Second, 
in a more realistic model, one must consider competition 
with DSC, which also could change features of the tran- 
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sition to IC state. A better understanding of IC DDW 
order parameter is necessary to better understand the 
the situation away from half-filling. 

III. C VS. IC ORDER 

A. C Order parameter 

Given that the susceptibility (|5J points towards an IC 
instability, the next issue is what IC DDW order would 
look like. Before addressing this question, let us first 
summarize some features of the C-casei& In the C case, 
a general singlet particle-hole condensate can be written 
as (neglecting the k = k' contribution): 

(C' t c Q , k )-*<f/(k)^'^, k+Q . (3) 

Here, $ is the magnitude of the order parameter. The 
wavevector 2Q is a reciprocal lattice vector. The spin 
index a included above will be omitted in the future for 
simplicity unless clarity requires otherwise. In J2J, the 
form factor /(k), which can transform non-trivially un- 
der the point group of the 2-D lattice, is related to the 
angular momentum of the particle-hole condensate. For 
DDW order, /(k) = cosfc^ — cosfc y which corresponds 
to I = 2 and d x 2_ y 2 symmetry. A p^-density-wave with 
Q = (0, 7r) has /(k) = sinfe^, which corresponds to / = 1. 
In real space, (0 becomes 

(4,c r }=V(r,r')$cos(Q-r) (4) 

where the vertices V of the orders discussed above are 
given by 

zV (r, r ) ^r',r+xa ^r'.r— xa ^r',r+ya ^r',r— ya 

-zV^ DW (r, r') = «5 r ,, r+xa - <5 r ,, r _ xa . 

(5) 

Note that the factors of i in V DDW and y p * DW signify 
that the corresponding phases break time-reversal, e.g. 
current flows along the bonds. Positive/negative signs in 
represent current going in/out of the vertex. Trans- 
lation invariance, by default, is broken by any density 
wave order since there is a preferred vector Q, but in the 
C case, invariance under a subset of the lattice group is 
retained. DDW order, for example, is invariant under 
translation by linear combinations of even reciprocal lat- 
tice vectors along each bond direction. The current along 
the bonds is given by: 

ir^v+sa = it((4c T+Sa ) - (ct +3a C T )) (6) 

Here s is along either the x- or y-direction. Because of 
the equivalence of Q and — Q in the C case, the cur- 
rent is simply j T <-*r+sa ~ cos(Q • sa) ~ (cjc r +£ a ). In the 
DDW case this results in an alternating plaquette cur- 
rent/checkerboard magnetic field patterns^ 



B. IC Order Parameter 

To obtain an IC version of DDW order, one first notes 
that DDW is chiral and hence breaks the Ising symme- 
try which reverses all the clockwise plaquette currents 
into counter-clockwise and vice versa. Therefore one can 
construct anti-phase domain walls as in Fig[5] A sim- 
ple caricature of a bond-oriented domain wall is given in 
FigEK- Its real space representation is 

(clc r ) = V DDW (r,r')$cosQ-rx(29(r :c )-l) + V IC (7) 

where Q(x) is the step function, so 

^ , x . / (2n + l)wx\ , . 
2e(a?)-l = 5^A2n+isin^ Tf^j ( 8 ) 

n * ' 

with A2 n +i = 4/7r(2n + 1). At the midpoint of the do- 
main wall, the currents in the x-direction vanish while the 
currents in the y-direction are halved. The diagonal do- 
main wall depicted in figure |2h has a similar representa- 
tion. Note that because of the halving of the y-direction 
currents at the midpoint of the bond-domain wall in J7J , 
we need to include the term: 

V IC = V R (I . r , )( ^ o + yL(r; (9) 

where the new vertices V L,R are given by: 

^ ^r',r-xa n^r',r+ya 77^r',r— ya 

l l (10) 

^r'.r+xa ^^r',r+ya ~^v' ,r— ya- 

Such vertices seem necessary for IC transitions of current 
density waves. One could imagine situations in which 
the transition to a disordered state is driven by thermal 
proliferation of such vertices, similar to earlier analysis 
of the 6-vertex model^. Although the analysis of such a 
transition is an interesting topic on its own, we will not 
pursue it in this paper. We are more interested in the 
energetics of the IC phase, so a long-wavelength mean- 
field treatment is sufficient. 

To analyze the C-IC transition, one then considers an 
array of alternating domain walls with average separa- 
tion L. Then, in the single harmonic approximation 
one retains only the n = term in JSJ and identifies 
q = tt / L with the incommensurability wavevector. The 
new short-range vertex contributions i|l(J|) are also ne- 
glected in the lowest harmonic approximation. The ne- 
glect of higher harmonics should not change the transi- 
tion temperature^ or the value of the critical doping for 
the onset of IC order since $ is small and the corrections 
are of higher order in <F An incommensurate instability 
in a clean system occurs when L becomes finite. (In a 
disordered system, it occurs when when L < Lg, where 
Ld is a characteristic length scale for the disorder.) 

Within the single harmonic approximation, an 
array of anti-phase domain walls with separation 
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q| = | Q — (7r, 7r) I = tt/L, looks simpler in momentum 
space: 

t * 

(Ck-Ck) = i-/(k)(<5 k /, k+Q + <5k/,k-Q). (11) 

Note that in the limit |q| — > 0, the above reduces to 
©■ One could have guessed the result (fTT|> by writing 
2<5k',k+Q = <5k',k+Q +<5k',k-Q in the commensurate case, 
and then let Q go incommensurate. This suggests other 
ways to get IC order. In particular the oddness of DDW 
order with respect to rotations by n/2 and transpositions 
about, say, the (n, 7r) direction lead to the following ex- 
ample of a DDW checkerboard pattern: 

(4'Ck) = J*/(k)(a*, k+Q + 5 k ,, k _ Q ) - {Q -> o(Q)} 

where the operation O(Q) can be either a transposition, 
or a rotaion by tt/2. The checkerboard pattern in 112|) 
corresponds to simply superimposing two domain walls 
rotated by tt/2 with respect to each other. 

In a similar spirit, one could also investigate the possi- 
bility of non-topological domain walls 30 , e.g. IC domains 
as the above, superimposed with an uniform order, so the 
overall order parameter does not change chirality. How- 
ever, as we will see in Sec. IIVI much of the physics of 
the C-IC transition will already be transparent from the 
simpler choices (|llf) and <|12fl of order parameter. 

Finally, note that within the single harmonic approxi- 
mation, and l|12fl are current-conserving only to low- 
est order in q. For the purpose of clarifying the energet- 
ics of the IC order close to the C-IC transition, or in the 
cases where the magnitude of the incommensurability is 
small, this lowest order approximation is good enough. 
Our numerical results, Sec. IIVI a posteriori show that 
we are in such regime. 

C. Microcopic Models 

We will consider the above order parameters in the 
context of the following Hamiltonian 26 : 

+ u ^2 n n n ji + v ^2 n i n j ( 13 ) 

In this formula, tij is hopping with — t for nearest 
neighbors, £y = t' for next nearest neighbors and fy = 
otherwise. t c is a correlated hopping term which simul- 
taneously hops an electron from site j to site i and hops 
an electron from i' into the vacated site j. The on-site 
and nearest-neighbor repulsions are, respectively, U and 
V. The indices i,j signify lattice sites and a the spin. 

We will consider this Hamiltonian at the mean-field 
level. At this level, the energetics of DDW order is the 




FIG. 2: (a) Bond-oriented DDW domain wall, (b) Diagonal 
DDW domain wall, (c) Current vertex of DDW. There is no 
current on dashed lines lines. Arrowless currents are regions 
of increased kinetic energy. Thick arrows represent one (ar- 
bitrary) unit of current. Thin lines represent half a unit of 
current. 



same as for a 'BCS-reduced' Hamiltonian for DDW order 
of the form: 

Hddw = -9ddw J /(k)/(k')cj c+Q(T c k(T c ] j c , (T ,c k ' + Q (T /, 

k,k' 

(14) 

with 5ddw = 24< c + 8V. A similar reduced Hamiltonian 
for d-wave superconductivity has gusc = I2t c — 8V. Note 
that the processes responsible for both DSC and DDW 
are essentially kinetic. 

A Hamiltonian of the form (|13|l has one important 
virtue: it has a d-wave superconducting ground state 
over a range of dopings (and an antiferromagnetic ground 
state at half- filling) . Since this is the sine qua non for 
any description of the cuprates, we believe that this is 
a good starting point. One may object that correlated 
hopping terms have to be put in by hand in a micro- 
scopic Hamiltonian. However, they naturally arise even 
from the one-band Hubbard model away from half-filling, 
where their contribution can be computed int/U pertur- 
bation theory: 

tc * ^ ( 15 ) 

At dopings x = 0.1, for t/U ~ 0.2, we would have 
ffDDW = 2grjsc ~ O.lt = J/2. Thus, not far away from 
half-filling, correlated hopping terms are certainly not 
negligible. 
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D. Mean Field Theory 

Let us generalize the mean field description of C DDW 
to the IC case. The path we will take is similar to that 
in Sec. IIII Bl where we wrote down the commensurate 
order parameter without making use of the Q — > — Q 
symmetry, and as a result the generalization to incom- 
mensurate Q was straightforward. We thus obtain a 
functional J-(/i, Ai,q), equation (|2UII . whose minimiza- 
tion yields self-consistently both the order parameters Aj 
and the deviation from commensurability q = Q — (jr, it). 
The only subtlety, as we discuss at the end of this section, 
is that although T reduces to the mean-field free energy 
in the C phase, in the IC phase it is not the free energy 
even within the single harmonic approximation. Instead, 
minimization of T is equivalent to a variational search of 
the ground state. 

Our derivations generalize the treatment of Nayak and 
Pivovaro\s2& of competition of C DDW with DSC and 
possibly AF. The C DDW mean field Hamiltonian we 
start with is: 



Hddw = J {^-^cLck+Qcr + Q 



keRBZ 



-Q} + {h.c.} 
(16) 

where W k = ^/(k). 

Note that in the Q — > — Q symmetry for in- 
commensurate Q is not equivalent to the Hermiticity of 
the Hamiltonian due to the k-dependance of Wk (e.g. 
(4ck +Q ) = zW k ? (c^Ck-q)* = -zT^k-q for IC Q). 
This is unlike the case of AF order, where the interac- 
tion is simply 

Waf = J 0(4 T c k+QT - c^Ck+Qx) +h.c. (17) 



keRBZ 



and we have (cj_T. tT<T 



kcr ' 3 Ck +Q(J <) = 4>= <cL T 3 CT<T ' c k-Q<T'}*- 

In order to check with earlier results on the C-IC tran- 
sition of AF order in the weak coupling Hubbard model^ 
(see Sec. |Hj|, we will also include a term like (|17|) in our 
full reduced Hamiltonian. Of course, in the strong cou- 
pling, U > t, limit, one is not justified in treating AF as 
a Fermi surface instability, especially close to half filling. 



At— J description is more appropriate in this case. How- 
ever, as mentioned in Sec. 1111 CI both DDW and DSC 
are essentially kinetic in origin within the extended Hub- 
bard model, so even in the strong onsite repulsion limit, 
there is a good justification to use mean field theory to 
describe their competition away from half filling. 

In the C case, the standard DSC term in our reduced 
Hamiltonian is: 



Hdsc = J Afc4t c -ki + h.c. 



(18) 



keRBZ 



where Ak = ^-/(k). Using a Nambu basis ^k = 

( c kt ' c iL+qt ' c — k i ' c — k— Qi. ) j the full Hamiltonian H = 
Wkin+^AF+WDDW+^DSC can be rewritten as H—fiN = 
Jrbz *k ' A • * k where: 



'£k~M «Gk 

^ ^ I c.c. e k +Q - /' 



c.c. 





c.c. 



A k 

A_ k -Q 

-£k + t 1 -iG-k-Q + <; 

c.c. -e k +Q + /i 



(19) 

with Gk = (Wk — W / k+Q)/2. The action derived from the 
above Hamiltonian yields, upon integrating the Fermions 
out, the free energy in the ordered state: 

^q = ^quad+Xl / { ( S£ k " - l ^ 2 COSh(^ ) } 

(20) 



^keRBZ 

with ±E S , s = ±1 the energy eigenvalues of l|19f) . and 



quad 



A 2 



5S DW ffDDW 3DSC 



(21) 



by: 



If we assume Q — > Q symmetry, E s would be given 



E a (k) = 



A k 



[M-e+(k) 



d(k)] 2 (22) 



where e±(k) = (e k ± ek+q)/2 and $ 2 = W 2 + </> 2 . 

However, in the IC case, Q — > — Q need not hold, in 
which case the above would generalize to: 



E.(k) = {f + (k)Al sc + (f i - e+ (k)) 2 + I- 2 + e 2 _(k) + syfifa - 6 + (k)) 2 ($ 2 + e 2 _(k)) + D^Al sc , $ 2 , £k , 
Z? q = [4/ k l>k + /£(k){(6k - m) 2 - (e k+Q - tf}]A 2 mc + £(k)A* S0 . 

I 



The effect of a nonzero incommensurability q can be 
seen not only in the appearance of the terms /±(k) = 

(/ 2 (k) ± / 2 (k + Q))/2 and / k = (/(k) + /(k + Q))/2, 
but also in the DSC-and doping-dependent function D q . 



Our final step before letting Q be incommensurate is 
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to pass the integration in (|20() from RBZ to BZ: 
■^-^ q uad= J f.(n,0,At) = ± J f B (n,0,Ai) (24) 

RBZ BZ 

Once we convert the RBZ integrals to BZ integrals as in 
Q24JI. we numerically minimize (see Sec. HVfl the resulting 
'free energy' with respect to both the order parameters 
Aj and the deviation from commensurability q. 




Doping x 



FIG. 3: Phase diagram of DDW competing with DSC. Both 
f/DDW = 0.05eV and <?dsc = 0.03eV are large (compare to 
Figs. [K] and [U) t = 0.5eV. The next nearest neighbor hop- 
ping, t' — — 0.05eV, is small compared to typical values fitted 
to ARPE$2i, but IC order was nevertheless robust away from 
half-filling. Inset: IC wavevector q at T — O.OfeV as a func- 
tion of doping. Wavevectors are in units of ir/a, where a is 
the lattice spacing. 

It must be emphasized that the expression to be mini- 
mized in Sec. HVl is not the true free energy of the system 
if Q is incommensurate. One way to see this is to note 
that the eigenvalues (|23[) are not invariant under k — > k 
if Q 7^ (jr, it) and therefore they cannot be interpreted as 
band energies. In a C-IC transition, the excitations split 
into a generally uncountably many bands even within the 
single harmonic approximation^. Therefore, the prob- 
lem of finding the mean-field free energy is notoriously 
difficult even in the single harmonic limit. The calcu- 
lation we are doing, on the other hand, is essentially a 
mean-field variational search for the ground state, and 
not for the excitations. A similar calculation was done 
by SchulaSi for the case of the C-IC transition of the 
weak-coupling Hubbard model close to half-filling. By 
continuity of the true free energy^, one expects that the 
resulting variational estimates for q = |q| would be accu- 
rate only close to the C-IC transition, where q is small. 



IV. NUMERICAL RESULTS 



bond DDW: 



diagonal DDW: V////////////A 




Doping x 



FIG. 4: Again jddw = 0.04eV, p D sc = 0.02eV are 'large' 
(note temperature scale) similar to Fig. [3] with t = 0.5eV. 
However this time we taka a more realistic t' = — 0.12eV. 
Again IC order is robust at dopings relevant for the pseudo- 
gap. Inset: IC wavevector q at T = O.OleV as a function of 
doping. See Fig. 1101 for a more detailed view of q (and also 
Fig. for the case of %' = -0.06eV ). 

To obtain the phase diagrams for the models in i|2U[l . 
for various fixed dopings x and temperatures T, we min- 
imize T = J-(Ai,q) with respect to both the order pa- 
rameters A, (as was done previously in the C-case2&) 
and the IC wavevector magnitude q. Each choice of IC 
order parameter (e.g. bond or diagonal) corresponds to 
a separate minimization procedure. To distinguish which 
IC order dominates at a particular point in (x, T) space, 
we compared the free energies of the orders considered, 
and chose the IC phase of minimal energy. Note that, 
as given in H2U|) . the free energy is a function of fx, so 
to obtain the dependance on the doping x, we first per- 
form a Legendre transformation and minimize instead 
J-{x) — J-(fi(x)) + n(x)N where the chemical potential 
is obtained from 1 — x = df /dfj, with / = T /Nq. 

Our approach is slightly different than that previously 
used in the C-case2&. There, the competition between 
AF, DDW and DSC was analyzed, and some sensitivity 
on the details of the microscopies was found. Phase di- 
agrams resembling that of the cuprates were obtained. 
The general trend was that any phase diagram which 
included d-wave superconductivity also had DDW order 
competing with it. By changing the bare coupling con- 
stants, rather different phase diagrams could be obtained 
and these generally did not have superconductivity. 

In our treatment here, we do not consider the compe- 
tition with AF order, which is not a Fermi surface insta- 
bility anyway in the context of the cuprate superconduc- 
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bond IC DDW: 




0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
Doping x 




0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
Doping x 



FIG. 5: t = 0.5eV, gnnw = 0.025eV, gusc = O.OleV. The 
next nearest neighbor hopping, t' = — 0.12eV. Note the ap- 
pearance of DSC proximate to 'magic' dopings. Inset: the IC 
wavevector q(x) at T = O.OOleV. See Fig llll for comparison 
between q(x) and n(_x). 

tors. In J20J, we still include AF, but only in order to 
check with previous results^! of IC AF order away from 
half filling in the weak coupling gAF < t limit. On the 
other hand, the AF order appearing at half filling in the 
cuprates is a strong coupling phenomenon. A reduced 
Hamiltonian (|17f) is the wrong starting point at half fill- 
ing in the limit U > t. 

Therefore, we only consider competition between 
DDW and DSC and look for IC order in this case. All 
the phase diagrams constructed this way are robust in 
the sense that DDW and DSC are found to coexist at 
mean field level for almost all reasonable values of the 
coupling constants. 

Sample phase diagrams are shown in Figs. In all of 
them, IC order is present at dopings 0.15 < x < 0.2 where 
DDW and DSC compete. In the first two of these figures, 
we vary t' at relatively high <7ddw an d <?dsc > while in the 
last two we explore the phase diagram for more realis- 
tic DDW and DSC couplings (compare the temperature 
scales to those in the cuprates). Note the complexity of 
Figs. 03 Figs. [BJcompared to Figs. Figs. 0] At realistic 
couplings, the IC wavevector profile is complicated, and 
furthermore DSC order appears in more than one doping 
region within the DDW dome. In fact, the appearance 
of incommensurability seems correlated with the appear- 
ance of DSC. A more systematic analysis to be discussed 
below reveals that generally the types of C-IC transi- 
tions that were produced in our phase diagrams fall into 
two categories: (a) those which correspond to motion of 
the bare Fermi surface away from commensurate nesting 
(shaded regions); (b) those which appear near certain 
'magic' filling fractions. Case (a) is similar to the famil- 



FIG. 6: Same as Fig. except jddw = 0.03eV. There is a 
distinct IC response (inset) at x = 1/8 and x = 1/16. See 
Fig. I12l for comparison between fi(x) and q(x). 




1 1 1 1 1 1 1 1 1 1 1 

0.00 0.01 0.02 0.03 0.04 0.05 

Temperature T/2t 

FIG. 7: Dependance of incommensurability on temperature 
for regions in the phase diagram in Fig. where DDW and 
DSC compete. Note that nonzero q is found above as well as 
below the DSC critical temperature. 

iar C-IC mechanism discussed by SchulaSI, while case (b) 
appears to be a consequence of negative effective stiffness 
at special fractions, the exact nature of which at present 
is not clear. In both cases q is independent of the onset 
of DSC, Fig. but because incommensurability relieves 
the DDW-DSC competition, the onset of DSC correlates 
to the presence of IC DDW. 

Let us first discuss (a). It is well knownSi that in the 
weak-coupling Hubbard model at t' = 0, the AF insta- 
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bility at half-filling is driven IC for any non-zero doping 
(see Sec. EJ)- The source of the C-IC transition can be 
traced to the approximate nesting of the bare Fermi sur- 
face (FS) at IC wavevectors away from half-filling. Close 
to half-filling, this picture would produce a linear devi- 
ation from commensurability, q ~ |/i — /i c ~ x — x c , 
Fig. [H] The incommensurability (deviation from (n,ir)) 
is horizontal or vertical, not diagonal, as we would expect 
because the bare susceptibility is more divergent due to 
a greater overlap between nested portions of the Fermi 
surface^!. 

A similar scenario happens in the case of DDW or- 
der. However, two notable differences from the AF case 
are present. First, DDW is less sensitive to t'/t, since 
it is driven mainly by antinodal regions of high density 
of states (d.o.s.). Therefore, while there is no AF insta- 
bility without nesting (at half-filling AF Fermi surface 
instability is only possible if qaf > 0(t')), DDW is still 
possible in this regime. This suggests that whenever the 
antinodal regions of the bare FS evolve away from the 
anti-nodal regons (which happen to also be regions with 
a high density of states as a result of the presence of van 
Hove singularities), a C-IC transition is likely to occur. 




r 1 i ' i ■ i 1 r — ■ — i — 1 — i — ■ — i — ■ — i — 1 — i — 

0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 



Doping x 



FIG. 8: AF at perfect nesting, t' = 0. The results of SchulaSi 
were recovered. Here, we used <?af = 2(7 = 0.05eV with t = 
0.5eV. Bond-oriented IC order was found to be energetically 
favorable, as a direct computation-- in the weak repulsion 
regime would yield. Inset: As expected, the IC wavevector at 
T — is linear in x. 

The second difference from the AF case is that DDW 
has nodes. This implies that the chemical potential in the 
ordered state need not be higher than the maximal DDW 
gap. As a result, DDW order can stay commensurate for 
wider ranges of doping. This is indeed what we find. For 
example, for all t'/t ~ 0, we obtain x^ ic < x™^ near 
zero temperature. In general, we find that the doping val- 
ues x c -i C for which the above type of C-IC DDW is possi- 



ble are unrealistically high in the context of the cuprates, 
if one considers more realistic bare FS, \t' /t\ > 0.25. Fur- 
thermore, in the IC DDW case, there need not be a simple 
relation q ~ |/x — fi c \ ~ x — x c between the doping and 
doping and deviation from commensurability. 

Let us now discuss the C-IC transitions that are 
present in our phase diagrams proximate to 'magic' fill- 
ing fractions x ~ 1/2". Although our understanding as 
to why they arise is incomplete, we would like to point 
out several interesting features about our results. First, 
to our knowledge, there is nothing special about the bare 
FS close to these filling fractions, likewise for the FS of 
the commensurate ordered state. The C-IC transition 
is not driven by nesting, since upon variation of t' we 
see similar signatures proximate to the 'magic' dopings. 
Second, the magnitude of the IC wavevector is generally 
small (corresponding to periodicity of more than 30 lat- 
tice spacings) and shows no clear dependance on n, Figs. 
151 4121 Third, we checked that similar behavior is present 
deep into the AF ordered stated, Fig. ED Unlike the 
DDW case, however, the IC lattice spacing in the AF 
case depends on n: q follows an increasing staircase pat- 
tern at 'magic' dopings. Note that our results for AF 
order seem to contradict Schulz's results^ which predict 
a linear dependence of q(x). 

The naive contradiction is easily resolved by realizing 
that Schulz's computational of q(x) relies on the Stoncr 
criterion, which is strictly valid at T c . Indeed, while deep 
into the ordered state the IC wavevector is stepwise in- 
creasing at 'magic' fractions, at higher temperatures we 
regain the linearity of q(x — x c ), Fig. Moreover, 
the fact that both the bare and the gapped FS in the C 
case are generally not nested at the the 'magic' fractions 
suggests that the transition is not driven by a divergent 
susceptibility, but is first order. Indeed, at low temper- 
atures, the chemical potential near the special dopings 
is non-monotonic for both the AF and the DDW+DSC 
diagrams. Similar behavior of fJ,{x) was found by Nayak 
and Pivovarov in the context of C DDW competing with 
DSG2&. 

A simple view of the energetics is given by the following 
Landau-Ginsburg free energy of two nodeless competing 
orders (the generalization to order parameters with nodes 
is straightforward): 

T= J d 2 x[p^\ 2 + Px {V X f + K{V 2 X f 

-m^\ 2 -m xX 2 + X\^\ 2 X 2 ) (25) 

In this free energy, \ plays the role of the DDW order 
parameter and tp the role of superconductivity. When 
p x becomes negative, \ orders at non-zero wavevector 
q = \J—p x /K. This occurs at the dopings indicated in 
Figures and EJ When ?p orders, this further favors in- 
commensurability if A > since the competition between 
the two orders is partially alleviated. This alleviation is 
a possible explanation of why the onset of DSC at low T 
is correlated to the presence of IC-DDW. 
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0.00 0.04 0.08 0.12 0.16 0.20 
Doping X 




0.00 0.05 0.10 0.15 0.20 0.25 
Doping x 



FIG. 9: Comparison between q(x) and ^(x) (inset) at T = 
0.01 for jddw = 0.04eV, g D sc = 0.02eV. Note that \t'\ = 
0.06eV which is smaller than the value taken in Fig. 0] 



FIG. 11: Comparison of /i(x) and q(x) (inset) at T — O.OOleV 
corresponding to the parameters in Fig. |K| Although here the 
IC response is not at 'magic' fractions, note that in Fig. 
there is DSC both at x — 1/8 and x = 1/16. 




FIG. 10: Same as Fig. EJexcept this time t' = — 0.12eV as in 
Fig. H 

The same phenomenon can be viewed from the oppo- 
site perspective: when DDW order becomes incommen- 
surate, there are regions in which it is suppressed (as 
in the domain walls of Fig. 2). Superconductivity can 
occur more easily there. Hence, incommensurability al- 
leviates the competition between DSC and DDW order 
at the mean-field level, thereby enhancing superconduc- 
tivity. Of course, IC order appears to suppress supercon- 
ductivity in the cuprates, but this could be due to effects 
beyond a simple mean-field theory. 

When the superconducting order parameter becomes 



FIG. 12: Same as Fig. 1111 (and at the same temperature) 
except this time there is a clear C-IC transition at 'magic' 
fractions (inset). At the same time /u(a?) is non-monotonic at 
these dopings. The coupling constants jddw = 0.03eV and 
ffDSC = O.OleV are the same as in Fig. 



strong enough, it may become energetically favorable for 
phase separation to occur. We can guess when this occurs 
from the usual Maxwell construction whenever the mean 
field chemical potential [J>(x) does not vary monotonically. 
This appears to occur at a few isolated 'magic' doping 
values. However, in a real system, proximity to a second 
order critical poinlsS, Coulomb repulsion^, or disorder 
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0.00 0.02 0.04 0.06 0.08 
Doping X 



FIG. 13: Example of the low temperature (T = 0.001 eV) 
staircase IC spectrum of AF around 'magic' filling fractions. 
The parameters we have taken here are (?af = — t' = 0.04 eV, 
t = 0.5 eV. Inset: fi(x) is non-monotonic at both jumps of 
q(x), but is smoothly decreasing whenever q is pinned. 



could drive the system in a state of mesoscopic phase 
separation. Moreover, the lattice pins any IC density 
wave order, so once the system is IC in order to minimize 
the free energy, the deviation from commensurability q 
would be pinned at a finite value. Finally it is important 
to remember that disorder could have important effects. 
It would not only disorder the array of domain walls, but 
could lead to droplets of one phase in regions of another. 

We stress that while the general trend towards IC or- 
der may be clear, many of the quantitative aspects are 
not clear, in particular the precise copings at which DDW 
order becomes incommensurate, it is not clear to us why 
such C-IC transitions happen near 'magic' filling frac- 
tions. The fact that similar behavior occurs in the case 
of AF order suggests that the tendency of phase sepa- 
ration at 'magic' filling fractions is a general trend for a 
variety of density wave order. Because the AF only case 
was explored upon taking the gDDW,DSC — * limit in 
the full free energy, it is not clear yet whether the ten- 
dency of phase separation at 'magic' filling fractions is 
purely a feature of density waves competing with other 
order parameters, or of the lattice. It might be worth 
reconsidering the AF case in the presence of other order 
parameters (such as DDW, DSC). What makes our re- 
sults surprising is that unlike in the frustrated phase sep- 
aration picture for stripe formation, there is no Coulomb 
repulsion present. Hence, there is no obvious connec- 
tion between phase separation and incommensurability. 
(Please note that while IC DDW order can induce charge 
density wave order at twice the wavevector, we have not 
included this effect, which we expect to be small, in our 
calculation.) Therefore, while Coulomb repulsion may 




O02 O04 O06 

Doping X 



FIG. 14: At higher temperatures q(x) regains its linearity 
around predicted by Schula^. Here T = 0.035eV and all 
other parameters are as in Fig. 1131 Inset: /x(x) is monotonic 
throughout, suggesting that at higher temperatures, the C-IC 
transition is driven purely by IC nesting of the FS away from 
half-filling. 

play a role in determining the CDW spacing at x = 1/8 
in LaCuO-i, in our calculation, the phase separation at 
'magic' fractions can only arise as interplay between the 
lattice, competing order physics and energetics. What 
exactly is the combination of the three deserves some fu- 
ture research^ 7 .. 



V. EXPERIMENTAL SIGNATURES 

According to the discussion following H25|) it seems 
likely that within the DSC state a competing DDW or- 
der parameter is IC as it develops. Hence, we briefly 
describe the experimental signatures that would distin- 
guish between C and IC DDW order. Originally, C DDW 
was also named 'hidden order' in the context of the pseu- 
dogap, because most of the experimental signatures as- 
sociated with it are indirect. The magnetic field, created 
by the alternating currents around neighboring plaque- 
ttes, is too small to be detected by neutron scattering 
(B ~ 1 — 30 G for Addw ~ 30 meV). At the same time, 
because of the symmetry of the C DDW magnetic field, 
there is no signature in NMR. The charge of the C DDW 
state is uniform so STM would not show signatures, un- 
less they are associated with interference of nodal quasi- 
particles due to scattering off impurities. It has been 
argued recently 3 ^ that such interference disperses weakly 
enough to account for the pseudogap d.o.s. incommensu- 
rate modulations observed by STM. For a further discus- 
sion of the signatures of C DDW, the reader is referred 
to the literaturo 7 . 10 i 13 i 17 i 18 . 19 i 20 i 21 . 
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An IC DDW is more easily detectable directly. The lo- 
cal electronic density of states would be modulated due 
to the IC order parameter, and hence could be observable 
by STM. Secondly, the IC modulations of the staggered 
currents would produce NMR line splitting of both the 
Cu and O atoms when the applied magnetic fields are 
perpendicular to the a-b planes. This is in sharp contrast 
with IC SDW, where there is an NMR line-splitting due 
to spatially modulated local magnetization, which is di- 
rected along the a-b plane. 

Let us estimate the modulated magnetic field strength 
arising from IC DDW (the reader is referred to2£ where 
similar considerations are applied to the case of orbital 
antiferromagnetism proposed to exist in U RU2S12) ■ Sup- 
pose we have an array of domain walls, Fig. along the 
y-axis. Deep between the domains the field is zero as in 
the case of C DDW. Along a domain at x = 0, where 
the field is maximal, again most of the fields cancel ex- 
cept those created by the vertical alternating currents 
of strength Iq at x = —a and those of strength Iq/2 at 
x = a. Equivalently this would be the field of a single 
vertical array of alternating currents of strength Iq/2 at 
x = —a: 

Bmax = ^(-1) 4 |(C0S(^ +1 ) - COs(^)) (26) 



where cos 6^ = [(i - l)a + yo]/[\/ ((* - l)a + Vo) 2 + Xq] 
and (xo,yo) = (a, a/2), (a, 0) for O, Cu atoms respec- 
tively. Clearly at yo — 0, Cu, atoms the field is zero by 
symmetry. For O atoms the sum can be performed yield- 
ing 4-7ri? max //!o = 0.14/o- In the dilute domain limit, this 
field is independent of the distance between the domains. 

For a sinusoidal current modulation, £? max would de- 
pend on the modulation length l q — 1/q. If I(y,x) = 
I(y)sin(qx) then a good estimate of B max is given by 
H26JI with Jo replaced by SIq = 2qaIo = 2Iq /n where n is 
the periodicity of the modulation (in units of the lattice 
spacing). Therefore, for large n we obtain 47r£? max //io = 
0.24/o/n. 

VI. CONCLUSION 

Incommensurate order is an important possibility 
when competing orders are considered. At a first-order 



transition the system is phase separated and hence nec- 
essarily inhomogeneous. If the competition between the 
orders is weak, coexistence is possible, but modulations of 
the competing orders are likely to be induced. The effects 
of the lattice, disorder, lomg-range Coulomb repulsion 33 
and proximity to criticality 32 is to generally stabilize in- 
commensurability on mesoscopic length scales. This is 
unlike the case of, say, the liquid to vapor phase tran- 
sition where the competition results in phase separation 
on macroscopic length scales. 

In the case of DDW and DSC, the competition is cer- 
tainly weak within the framework of the extended Hub- 
bard models we have considered. Because the processes 
stabilizing both orders are of similar nature and are at 
the same time smaller than all of onsite repulsion, ex- 
change and hopping, it does not cost much energy to 
convert DDW order into DSC and vice versa. A sim- 
ple resolution of the competition is for DDW order to 
become incommensurate. We find that this occurs for 
some doping levels, but the incommensuration tends to 
be relatively small. Such IC DDW could, in principle, be 
observable at low temperatures by NMR since it creates 
an inhomogeneous magnetic field at O atoms perpendic- 
ular to the ab-plane. Most of NMR experiments test for 
much stronger magnetic fields parallel to the ab-plane on 
the Cu atoms, and at higher temperatures such IC DDW 
phase would be smeared by thermal fluctuations, which 
may exlpain why such IC order, if it exists, it hasn't been 
seen by experiment. 
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